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AN HST/COS SEARCH FOR WARM-HOT BARYONS IN THE MRK421 SIGHTLINE 
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ABSTRACT 

Thermally-broadened Lya absorbers (BLAs) offer an alternate method to using highly-ionized metal 
absorbers (O VI, O VII, etc.) to probe the warm-hot intergalactic medium (WHIM, T — 10 5 — 
10 7 K). Until now, WHIM surveys via BLAs have been no less ambiguous than those via far-UV 
and X-ray metal-ion probes. Detecting these weak, broad features requires background sources with 
a well-characterized far-UV continuum and data of very high quality. However, a recent HST/COS 
observation of the z = 0.03 blazar Mrk421 allows us to perform a metal-independent search for 
WHIM gas with unprecedented precision. The data have high signal-to-noise (S/N « 50 per ~ 
20 km s _1 resolution element) and the smooth, power-law blazar spectrum allows a fully-parametric 
continuum model. We analyze the Mrk421 sight li ne for BLA absorbers, particularly for counterparts 
to the proposed O VII WHIM systems reported bv lNicastro et all (j2005al fb1) based on GTiancira/LETG 
observations. We derive the Lya profiles predicted by the X-ray observations. The signal-to-noise 
ratio of the COS data is high (S/N 25 per pixel), but much higher S/N can be obtained by binning 
the data to widths characteristic of the expected BLA profiles. With this technique, we are sensitive 
to WHIM gas over a large (Njj, T) parameter range in the Mrk421 sight line. We rule out the claimed 
Nicastr o et al.l O VII detections at their nominal temperatures (T ~ 1 - 2 x 10 6 K) and metallicities 
(Z = 0.1 Zq) at > 2(j level. However, WHIM gas at higher temperatures and/or higher metallicities 
is consistent with our COS non-detections. 

Subject headings: intergalactic medium, quasars: absorption lines, cosmology: observations, BL Lac- 
ertae objects: individual (Mrk421) 
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1. INTRODUCTION 

The Warm-Hot Intergalactic Medium (WHIM) is ex- 
pected to make up 25-50% of the baryons in the 
low-re dshift Universe based on cosmological si mula- 
tions (iCen fe Ostrikerl H99l iDave et all [20011 
ISmith et al.l 201lfh Thus far, only a small fraction o f 
this hot gas has been detected (see, e.g. lDanforthll2009l) . 
and the interpretation of these detections is still contro- 
versial. The difficulty in detecting this gas lies in its com- 
bination of high temperature (10 5-7 K) and low hydro- 
gen density (10~ 4 — 10 -6 cm -3 ). For example, only the 
subset of WHIM gas which is sufficiently metal-enriched 
can be detected through highly-ionized metal ions at UV 
and X-ray energies. The O VI AA1031, 1037 doublet has 
been detected in numerous extragalactic systems with 
the Far Ultraviolet Spectroscopi c Explorer (FUSE) and 
Hubb l e Space Telescope (HST) (IDanforth fc Shulll 12001 



l2?M iTriPD et all 120081: iThom fc Chenl 120081 ). It may 

trace very low density, photoionized ~ 10 4 K gas (e .g., 
lOppenheimer fc Dave 2008; Qpp enheimer et al. 2 01l[) or 
shock -heated, collisionallv-ionized WHIM (jSmith et all 
1201 lh . Until now, the searches for hotter WHIM (> 
10 6 K) have been restricted to soft X-ray absorption. 

Very highly ionized metal species such as O VII, O VIII, 
N VII and Ne IX with soft-X-ray resonance transitions 
offer a good option for finding diffuse gas at the hotter 
end of the WHIM temperature range. However, the pre- 
dicted low metalli city of this gas (Z ~ .1 Zq typical 
of O VI absorbers: IDanforth fc Shulll 120081 ) and the rel- 
ative insensitivity and low resolution of even our best, 



current X-ray spectrometers (i.e., the Chandra Low En- 
ergy Transmission Grating (LETG) and XMM-Newton 
Reflection Grating Spectrometer (RGS)) has generated 
controversy for most or all current proposed detections 
of hotter WHIM gas. 

The growing number of "broad Lya absorbers" (BLAs) 
detected in h igh si gnal-to-noise F UV s pectra 
dSembach efafl (2001 iRichter et all l200l 120061: 
iLehner et al.l 120071; IDanforth et al.l l2010al) may trace 
WHIM independent of metal enrichment, but this 
diagnostic suffers from a different set of observational 
uncertainties. Even at the cooler end of the WHIM 
temperature range, the fraction of hydrogen in a 
neutral state is tiny. Thermal broadening renders any 
Lya absorption shallow and broad, and the WHIM 
absorber may be blended with narrower absorption 
features arising in adjacent, cooler, photoionized gas. 
At temperatures above 10 6 K, spectra of extremely high 
quality, with well-defined continua, are required to char- 
acterize the WHIM through BLAs alone. Independent 
measurements of individual WHIM systems through 
complimentary techniques reduce the uncertainties 
inherent in single-method detections. 

Two of the earliest and most-cited X-ray WHIM 
detections wer e the two reported O VII absorbers 
(|Nicastro et al.l l2005allH ) in Chandra LETG observa- 
tions of the blazar Mrk421 (z em = 0.030). These 
observations were made while Mrk421 was in a flar- 
ing state and represent some of the highest signal-to- 
noise soft-X-ray data taken to date on an extragalactic 
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Fig. 1. — HST/COS spectrum of Mrk421 spanning the entire redshift pathlcngth of the Lya forest (1216A < A < 1252 A). The data in 
the top panel are binned by three pixels (~ 1/2 resolution element). The power-law continuum, Galactic damped Lya absorption, intrinsic 
Lya emission (z em = 0.030), and known narrow absorption features detailed in the text are modeled by the red curve offset from the flux 
by +5 X 10 — 14 erg cm -2 s _1 A -1 . The lower panel shows the residual data binned by seven pixels (approximately the COS resolution 
element) after the continuum and all absorption/emission features are divided out. The regions where narrow lines have been modeled out 
of the data are shown by orange bars. Blue diamonds mark the redshift locations and ±300 km s _1 uncertainties of the INicastro et al.l 
(2005a) claimed O VII features. 



sight line. Although these absorber s were not confirmed 
in observations by XMM-Newton flKaastra et all 120061 : 
iWilliams et all 120061: iRasmussen et al.ll2007t) . their pres- 
ence has neither been confirmed nor ruled out conclu- 
sively. While a few other X-ray WHIM absorption de- 
tections hav e been proposed (e.g., in the PKS 2155-304 
sight line bylFang et aLll2002D . a new tactic employed by 
IBuote et al. (2009) may have finally achieve d some suc- 
cess. Based on the results of simulations (e.g. JDave et al.l 
1999) that show high concentrations of WHIM gas 
aroun d large-s cale filamentary structure. IBuote et al.1 
(2009) (see also lFang et al.l^OlOtrZappacosta et al.ll2010| ) 
observed the bright blazar H 2356—309 behind the Sculp- 
tor Wall and made O VII detections at the redshift 
of that filament with both Chandra and XMM-Newton. 
Based on the Sculptor Wall detection, it seems worth- 
while to revisit the Mrk421 sight line, since one of the 
two claimed WHIM detections is at the redshift (z = 
0.027) of anothe r nearby, large-scale galaxy filament , 



the Great Wall (Ide Lapparent. Geller fc Huchraj 19861 : 
iPenton. Stocke. fc Shulll)2000l : IWilliams et al.ll20iot )~ 



In this paper, we focus on the potential for detecting 
WHIM in high-quality UV spectra from the Cosmic Ori- 
gins Spectrograph (COS) on HST. This analysis is spe- 
cific to the Mrk421 sight line, but easily generalized to 



sight lines for which high-quality far-UV spectra exist. 
The FUV spectra of blazars can be well-mod eled as a 
power law ([Stocke. Danforth. fc Perlmanll201ll ). and the 
high S/N and multiple ionic transitions allow foreground 
Galactic absorbers to be well-modeled. Since the data 
can be fitted with a fully parametric continuum model, 
this is an ideal test-bed for detecting WHIM gas in BLAs. 

In Section 2, we present the COS FUV spectrum 
of Mrk421, identify all spectral features present, both 
Galactic and extragalactic, and produce a normalized, 
line-free spectrum covering the entire Lya forest region. 
In Section 3, we model the BLA absorb ers implied by 
both the reported INicastro et al.l (|2005al fbT) (henceforth 
collectively N05) O VII measurements and for more gen- 
eralized WHIM absorbers. We then place limits on the 
temperature and metallicity of any WHIM gas present 
along the Mrk421 line of sight based on the observed 
COS data. In Section 4 we discuss the meaning and im- 
portance of these results, reiterate our main conclusions, 
and discuss applications to other sight-lines. 

2. OBSERVATIONS AND ANALYSIS 

COS far-UV observations of BLLac object Mrk421 
(z em = 0.030) were carried out in 2009, December, 
during the first three months of COS science obser- 
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vations as part of the COS Guaranteed Time Obser- 
vations (PID 11520, PI Green). Four exposures were 
made in the G130M (1135 A < A < 1480 A; totalling 
1.7 ksec) medium-resolution grating (R ss 18,000), 
each at a different central wavelengt h setting, to dither 
over known instrumen tal features (see [Green eTal]l20ll 
lOsterman et alJl201lD . 

The exposures were reduced with a custom flat field 
derived from the combination of several public COS 
calibration observations of white dwarfs. Four targets 
were used to create the G130M flat field; WD 0308-565, 
WD 0320-539, WD 0947+857, and WD 1057+719. For 
each calibration target, an iterative technique was used 
to independently measure the relatively smooth con- 
tinua from the inferred underlying flat-field and grid wire 
shadows using the 'xld' files output by CALCOS 2.13. 
The four independent one-dimensional estimates of the 
G130M flat field were merged in detector space with 
(S/N) 2 weighting. After applying the custom flat field 
to each exposure, we aligned the four Mrk421 exposures 
in wavelength space via a cross-correlation of regions 
around prominent ISM absorption features. Exposures 
were interpolated onto a common wavelength grid and 
combined via an exposure time weighted mean as de- 
scribed in lDanforth et al.l (|2010bl) . 

Our custom flat field is not perfect, but it consistently 
improved the continuum signal-to- noise in our Mrk421 
spectru m as compared t o the standard "wireflat" pro- 
cessing (jDanforth et al.l l2010bf ). For example, in the 
1410-1430 A region of the Mrk421 G130M spectrum, 
where no ISM or IGM absorption lines are expected, 
the signal-to-noise per seven-pixel resolution element im- 
proved from 27 to 36 with our custom flat fielding, as 
compared to our standard "wireflat" processing. A sec- 
ond line-free region from 1340-1369 A increased from 
S/N of 38 to 48 per resolution element. 

2.1. Spectral Modeling 

The Lya forest region of the Mrk421 sight-line is 
shown in the top panel of Figure 1. Since we are inter- 
ested in weak, broad absorption features in the data, an 
accurate and detailed model fit to the spectrum including 
all known features is crucial. Fortunately, as with most 
blazars, the dereddened UV continuum is easily defined 
by a power law of the form Fx = F 912 (A/912l)- Q \ We 
use the parameters Fq = 3.405 x 10~ 13 erg cm -2 s _1 A" 1 
and ot\ = 1.918 ± 0.020 to provide a baseline contin- 
uum fit to the data in the 1150 A-1300 A region. Note 
that these parameters are slightly different from those 
reported bv lStocke. Danforth. fc Perlmanl ()201lD . which 
were based on an earlier reduction of the data and a con- 
tinuum fit to the entire COS/FUV spectral region. The 
damped Galactic Lya absorber is modeled with the pa- 
rameters v — —28 ± 1 km s _1 , &hi = 55 + 2 km s _1 , and 
log N m (cm~ 2 ) = 20.042 + 0.002. This is consistent with 
the emission- weighted mean H I 21 cm profile in a ~ 0.6° 
beam toward Mrk421 in the Leiden - Argentine-Bonn H I 
survey (LAB. iKalberla et al.ll2005h . The measured col- 
umn of neutral hydrogen corresponds to a Galactic red- 
dening of E(B - V) = iV ffl /5.8 x 10 21 cm" 2 = 0.019. 
(A fit to the Galactic Ly/3 profile in the FUSE data 
yields a consistent value.) In a departure from the pure 
power-law source continuum typical of BLLac objects, 



Mrk421 shows a weak emission feature at 1251.9 A cor- 
responding to Lya emission at z em = 0.0298 from cir- 
cumnuclear gas in the blazar . We include the fit from 
iStocke. Danforth. fc Perlinanl (poTlh (FWHM = 1.23 
A, I = 1.27 x 10~ 14 erg cm~ 2 s _1 ) in our spectral model, 
shown offset from the data in the top panel of Figure 1. 

There are several narrow Galactic absorption lines 
which can be modeled out of the spectrum. The stronger 
blue line of the N V AA1238.82, 1242.80 doublet is well 
fitted with parameters v = —21 ± 6 km s _1 , b — 59 ± 8 
km s _1 , and log 2Vnv = 13.34+ 0.05, but the weaker 
red line of the doublet shows an absorption excess that 
is discussed below. The S II 1259.52, 1253.81, 1250.58 A 
triplet has a strong Galactic component (v — —8 + 2 
km s-\ b = 13 ± 4 km s" 1 , log N su = 15.05 ± 0.10), a 
weaker, blue-shifted component (v — —48 + 2 km s _1 , 
b = 14 ± 5 km s _1 , log iV s „ = 14.89 ± 0.05), and 
very weak red wing (v = 32 km s" 1 , b = 20 km s" 1 , 
log Nsu = 13.9 ± 0.1). These low-ion velocity compo- 
nents are consistent with the two H I emission compo- 
nents (|Wakker et alj|2003 IKalberla et aT]|2"005l ) and ab- 
sorption seen in o ther low-ionization species in COS data 
(|Yao et all 1201 11 ). The strong Si II A1260.42 line, while 
redward of the Lya forest region in Mrk421, can be fitted 
with two components (v — —21 ± 3 km s^ 1 , b — 28 ± 3 
km s _ \ log Nsai = 15.2 ± 0.1) and (v = -117 ± 5 
km s _1 , b = 14 ± 5 km s _1 , log N sill = 13.1 ± 0.1). 
The last identified Galactic absorber is the weak, blended 
Mg II AA1239.93, 1240.40 doublet. Assuming a 2:1 ra- 
tio of equivalent widths and similar velocity centroids 
and 6-values to the other Galactic lines, we model the 
system with v = —21 km s _1 , b = 59 km s _1 , and 
log iV M gn = 15.3. 

The sight line toward Mrk421 is relatively short (Az = 
0.030), and there is only a s ingle previously-reported 
(Shull. St ocke. fc Pentonl Il996h intervening Lya forest 
system (z = 0.01014, A = 1228.0 A). This moderate- 
strength Lya line is well-fitted by a Voigt profile with 
parameters (z = 0.01013, b = 19 ± 1 km s -1 , W\ = 
76 ± 4 mA, log N m = 13. 25 ± 0.02) convolved with the 
COS line spread function (jGhavamian et al.l l2009: Kriss 
1201 lh . No corresponding metal absorption (e.g., C IV, 
Si III, etc.) is seen at this redshift. 

A weak absorption feature (W\ = 25 ±4 mA) is seen at 
1246.4 A in the coadded spectrum. There are no species 
commonly seen in the ISM with absorption lines near 
1246 A, and the normalized spectrum in the 1243 — 1250 
A range is free of line-removal artifacts. The most likely 
identification of this feature is a weak Lya forest system 
with fit parameters cz — 7582 ± 8 km s _1 (z — 0.0253), 
b = 52 ± 12 km s" 1 , and log N m = 12.7 ± 0.1. 

The spectrum is divided by the model (offset curve 
in the top panel of Figure 1) to give a normalized spec- 
trum, with all identified absorption and emission features 
removed (Figure 1, lower panel). With the exception of 
a few narrow line-removal artifacts, the normalized spec- 
trum shows a remarkably uniform appearance over the 
range 1218 A < A < 1257 A. 

3. WHIM ABSORPTION TOWARD MRK 421 

Given a few input parameters, we can model the ex- 
pected appearance of a BLA. These parameters can be 
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TABLE 1 

Mrk421 WHIM Absorber Parameters 



Quantity 



0.011 



■■ 0.027 



Units 



Source a 



UV/X-ray Measurements 



cz 


3300 ± 300 


8090 ± 300 


km s 


-l 


N05 


No vn 


1.0 ±0.3 


0.7 ±0.3 


10 15 


cm z 


N05 


Novi 


< 1.5 


< 1.5 


10 13 


_9 

cm 


S05 


Ns vii 


0.8 ±0.4 


1.4 ±0.5 


10 15 


_9 

cm 


N05 


logT 


6.13±0.39 b 


6.15 ± -0.15 c 


K 




N05 



Implied Lya Absorber Parameters 





14S+H 


152ig 


km s 


-1 


log 7V m d 


12 5 +a6 


12.7 ±0.2 


z 0.l 


cm - 


w Lyoi d 






7 -l 


mA 


TO 


0.016 


0.025 






SL 


1.9 


2.7 


£7 





Savage et al. 2005 



a N05 = Nicastro et al. 2005a (Chandra); S05 
(FUSE) 

b Inferred from O VI/VII/VIII ratios 

c Inferred from NVI/VII ratios 

d Assuming Z = 0.1 Zq. See text for details. 



assumed for WHIM absorbers in general, or they can be 
derived from X-ray WHIM measurements and a few basic 
assumptions. For general single-phase WHIM gas with 
total hydrogen column density Nh at temperature T, the 
neutral hydrogen column density is Nm — fm(T) Njj 
and the thermal line width is 



b(T) = (40.6 km s" 1 )^ 



1/2 



(1) 



where T5 is temperature in units of 10 5 K. The neu- 
tral fraction as a function of temperature, /hi(T'), can 
be determined either in collisional ionization equilibrium 
(CIE) or with a more sophisticated sets of assumptions. 
In the following analysis, we as sume Z = 0.1 Z Pl gas coo l- 
ing isobarically as tabulated in lGnat fc Sternberg (|2007| ). 
At WHIM temperatures, the hydrogen neutral fraction is 
identical to that in CIE and the high-ion fractions differ 
only slightly from their CIE values. 

If we start from the reported WHIM parameters de- 
rived from the Chandra observations, we can derive the 
Nhi via the additional relationship, 



N m 



/hi(T) JVqvii 



(2) 



/ovn(T) (0/H) Q Z ■ 

Taking the reported INicastro et al.l X-ray detections at 
face value, the cz = 3300 km s _1 X-ray absorber is as- 
sumed to have T ranging from 5.5 x 10 5 K to 3.3 x 10 6 K, 
based on O VI/O VII/O VIII column density detec- 
tions and limits. The cz = 8090 km s _1 system is in- 
ferred to have a somewhat narrower temperature range, 
T = 1-2x1 6 K from N VI/N V II ratios and an 
O VI non-detection (jSavage et al]|2005D . These WHIM 
parameters predict Lya absorbers with column density 
log N m = 12 A - 13.3, b = 95 - 233 km s" 1 , and W x = 
13 — 100 mA at the z = 0.011 absorber; and log Nm w 
12.7, b = 130 - 180 km s~\ and W\ = 10 - 34 mA at 
z = 0.027 (Table 1). We assume (0/H) Q = 4.90 x 10~ 4 
(jAsplund et al.l I2009D and Z = 0.1 Z Q . Lower abun- 
dances of oxygen would scale these predictions upwards. 
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Fig. 2. — The detection limits for WHIM gas via broad Lya ab- 
sorption depend on the parameters of the absorber and the quality 
of the far-UV spectrum. In HST/COS data with (S/N) 1 » 25 
per pixel, idealized, unblended BLAs can be detected at 1 — 5 cr 
significance levels as shown by the dotted contours. Typical de- 
tection limits for several different column densities of O VII (blue, 
solid) and O VI (red, dashed) are shown for an assumed metallicity 
of Z = 0.1 Zq. Metal line diagnostics scale as Z^ 1 on the verti- 
cal axis. Blue diamonds sh ow the parameter-space location of the 
two Nicastro et al. (2005a) claimed WHIM detections with their 
uncertainties. 



3.1. Signal-to-Noise and Significance Levels 

The predicted BLAs are very weak (line-center optical 
depth To < 3%) and much broader than the ~ 18 km s _1 
instrumental resolution of COS. They will be difficult 
to detect in even high-quality HST/COS spectra and 
will require a good understanding of the signal-to-noise 
ratio and systematics in the data. The per-pixel scat- 
ter in the normalized, coadded data gives signal-to-noise 
per pixel (S/N)i = o~ x = 25 ± 1 over the Lya for- 
est wavelength range. In the Poissonian ideal, pixel- 
to-pixel noise is uncorrelated, and combining adjacent 
pixels results in lower noise (and hence higher S/N) by 
the square root of the number of pixels binned. We 
find that binning or smoothing the line-less, normal- 
ized data by x pixels results in an empirical relation- 
ship (S/N) x = (S/N)i x o.38±o.02 ; consistent witn the be- 
havior in other data sets over a wider wavelength range 
(|Keenev et al.l [20TlT ) . The fact that the index on x is 
less than 0.5 shows that the noise in adjacent pixels is 
correlated to some degree, i.e., fixed pattern noise exists 
in the data. 

The significance level of a detection depends on the ob- 
served equivalent width W\ and the S/N of the data. Fol- 
lowing the method of iKeenev et al.l (|2011l) (their equa- 
tion A6), we calculate the significance level, SL, of an 
absorption feature in data binned by x pixels as 



SL(x) = 



f c (x) W x (S/N) x _ f c (x) W x (S/N), x 



u.:is 



AA, 



x AAq 



(3) 

where f c (x) is the fraction of the total line profile en- 
compassed by x pixels. For lines much broader than the 
instrumental resolution. IKeenev et all find that SL(x) is 
maximized when the data are binned to approximately 
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Fig. 3. — Predicted broad Lya absorbers superimposed on the normalized, line-less COS spectrum of Mrk421. In the top panel, the 
COS data are binned to seven pixels (<~ 1 resolution element, ~ 17 km s — 1 ). Orange bars show regions where narrow absorption features 
have been modeled out of the spectrum. The blue curves offset below the data show the range of BLA absorption profiles pre dicted from 
the O VII detections (solid blue lines) and ±lcr uncertainties in T and JV(0 VII) (dashed blue lines) of Nicastro et al. (2005a) as listed in 
Table 1. Blue diamonds and horizontal error bars show the velocity centroids and ±300 km s _1 uncertainty of the Chandra detections. 
The bottom panel shows the improvement in (S/N) x obtained by binning the data by 50 pixels. Optimal binning for BLAs at T ~ 10 6 K 
is ~ 100 pixels. 



the FWHM of the line or i « 2Vhi2 6 and that this 
binning encompasses f c (x) = 76% of the total line area. 
Since each COS pixel in the region of the Lya forest is 
approximately 2.4 km s , this means 

SL(b) = 0.095 x W x , (4) 

with W\ in units of mA and 6 in km s — 1 . 

Applying Eq. [4] to observed lines, we confirm that the 
Lya forest line at 1228.0 A is highly significant (SL ~ 
30a). The weak absorption feature at 1246.4 A assumed 
to be Lya at z — 0.0253 is a ~ 5c detection in optimally- 
binned data. Similarly, the stronger (1238.82 A) line of 
the Galactic N V doublet is a ~ 9a detection and the 
blended Galactic Mg II doublet is significant at ~ 3a. 

3.2. WHIM Parameter Space Accessible by BLA 

Figure 2 shows the temperature and total column den- 
sity parameter space in which we might expect WHIM 
systems. HST/COS data with {S/N) 1 = 25 are sensitive 
to BLAs in the upper left half of the parameter space. 
Also shown are the regions in which WHIM could be 
detected via X-ray (O VII) and UV (O VI) metal-line di- 
agnostics at various column densities in Z = 0.1 Zq. The 
metal-line sensitivity limits scale on the vertical axis as 
Z^ 1 gas. 

Can we confirm or disprove the X-ray WHIM detec- 
tions based on the COS data? Figure 3 shows the COS 
data with identified absorption features removed (see 
Fig. 1) and then binned to 7 and 50 pixels in the upper 
and lower panels, respectively. The WHIM absorbers 
predicted from the X-ray measurements would appear 
at A - 1230 A and A ~ 1249 A (blue diamonds), both 
with a velocity uncertainty of ±300 km s _1 from N05. 
As predicted, broad features appear more significant in 
the lower panel where the pixel-to-pixel noise has been 



reduced and the binning width is closer to the charac- 
teristic line width. The N05 WHIM absorbers modeled 
via Eq. 3 are predicted to produce BLAs with doppler 
widths b « 150 km s~ x and equivalent widths of ~ 18 
and ~ 26 mA for the z = 0.011 and z — 0.028 features, 
respectively (blue curve in Figure 3). From Eq. HJ these 
would be significant (modulo the considerable range in 
predicted BLA properties; see Table 1) at the 1.9cr and 
2.7ct levels in (S/N)i = 25 HST/COS data binned to 
~ 100 pixels. 

It is clear from Figure 3 that no ~ 2 — 3a absorption 
features appear in the COS data consistent in either line 
profile or centroid with the predicted WHIM absorbers. 
The claimed Great Wall WHIM detection should appear 
at 2 — 3a significance at 1248.5 ± 1.2 A, but is not seen 
in the COS data. The 1246.4 A feature was discussed 
above and was removed from the residual spectrum in 
Figure 3. While relatively broad for a Lya feature (b = 
52 ± 12 km s" 1 ; T < 10 5 K), it is inconsistent with the 
iNicastro et all O VII detection at cz = 8090±300 km s" 1 
in both redshift and inferred temperature. The 1228.0 A 
Lya line is consistent in redshift with the claimed z — 
0.011 O VII system, but is far too narrow to arise in gas 
at T ~ 10 6 K. 

We therefore rule out WHIM absorption for the two 
systems wi th the redshift. tempe rature and metallicity 
reported bv lNicastro et al.l (|2005a| ) at a > 2a level. How- 
ever, if lower O VII column densities, higher temper- 
atures, or higher metallicities are assumed, the (Njj,T) 
parameter space locations of the claimed WHIM systems 
move toward the lower right in Figure 2 and would not 
be detectable in the current COS data. 

Two additional weak features become apparent in Fig- 
ure 3 in locations unrelated to the O VH-predicted BLAs. 
While both appear significant in the highly-binned lower 
panel of Figure 3, they are best fitted as single, narrow 
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components more typical of the photoionized Lya for- 
est. The first is a feature near 1225.8 A, a region uncon- 
taminatcd by narrow absorption lines. It is well-fitted 
as a ~ 4<7 narrow Lya system (b = 25 ± 11 km s _1 , 
W x = 13 ± 4 ml, log N m = 12.4 ± 0.1) at z = 0.0083. 
The second residual feature appears as an absorption ex- 
cess in the weaker Galactic N V line at 1242.5 A. Mea- 
surement of this feature is complicated by the overly- 
ing N V absorption, but if the absorption profile of the 
stronger N V line is applied to the weaker line, the resid- 
ual can be fitted as a ~ 5cr feature: b ~ 27 km s -1 , 
W x ~ 19 ml, log N m ~ 12.6, z = 0.0220 if interpreted 
as Lya. 

4. CONCLUSIONS AND IMPLICATIONS 

The INicastro etaT] (|2005al) claims of O VII absorbers 
toward Mrk421 were met with a great deal of skep- 
ticism in the X-ray c ommunity (Kaa stra et al.l 120061 : 
iRasmussen et all [2007). Even the more significant z = 
0.027 O VII system in the Great Wall (AWn ~ 7 x 
10 14 cm" 2 ) is only a ~ 2 mA feature in the Chan- 
dra LETG data and relies critically on the continuum 
model in the 21-23 A region. While Mrk421 has 
been observed for a total of > 200 ksec over many dif- 
ferent observing campaigns and features some of the 
highest S/N data for any extragalactic X-ray target, 
a more reliable O VII detection limit may require 
(> 5 mA). Note that the strong 
~ 2 ) O VII absorption at z — toward 



A^ovn > 3 x 10 15 cm" 2 
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cm 



(JVovn ~ 10 

Mrk421 is uncontroversial. Indeed, it is typical of z = 
absorption s een in numerous extragalactic sigh t lines 
dFang et al.l 120021: iFang. Sembach. fc Canizaresl 120031: 
McKernan. Yaaoob. fc Reynolds] 120041: INicastro et all 
20021 12005a): IWang et al.1 120051: IWilliams et al.1 120061: 



Buote et al 



20091) and is thought to arise in hot gas in 



the Galactic halo. 

In this paper we highlight the diagnostic power of 
thermally-broadened Lya absorbers in far-UV spectra as 
tracers of WHIM gas. Mrk421 was observed for less than 
half an hour (1.7 ksec or half an HST orbit) in the blue 
G130M mode of the Cosmic Origins Spectrograph, and 
yet it reaches WHIM detection levels comparable to those 
in > 200 ksec of X-ray observations with Chandra. Based 
on these high-quality data, we can rule out the existence 
of WHIM gas over a large region of (T,Njj) parame- 
ter space (see Figures 2, 3) along much of the Mrk421 
sight line. Specifically, we rule out the existence of broad 
Lya counterpart absorption to the claime d O VII WHIM 
absorb ers at z = 0.011 and z — 0.027 (Nicastro et al. 
l2005al[ bT) at ~ 2a level. We cannot, however, rule out the 
existence of the claimed WHIM absorbers if their metal- 
licities are higher than the assumed Z > 0.1 Zq, if the 
O VII column densities are at the lower limits of their re- 
ported ranges, or if the temperatures are higher than the 
~ 1 — 2 x 10 6 K inferred fr om ion ratios (jNicastro et al.l 
I2005al|bl : [Savage et alll20"05l ) . The Lya absorption m any 
of these systems would then be below the detection limit 
of the current COS data. 

The search for missing baryons has long centered on 



the gas in the WHIM phase via one of several tech- 
niques. Roughly one hundred \ow-z O VI detections 
have been reported in ~ 30 sight-lines (iDanforth fe Shulll 
[20081: iTripDeTalll^oM PThom fc Chenll2008n which give 
a statistically-significant sample of lower temperature, 
higher metallicity WHIM baryons (T < 10 6 K; Z > 
0.1 Zq; see red contours in Figure 2). However, while 
the presence of these absorbers is agreed upon in most 
cases, their interpretation as WHI M gas is still co ntrover- 
sial in roughly half the cases (see, Danforth 2009). X-ray 
metal lines such as O VII seem to present a more se- 
cure identification of WHIM absorption, but suffer from 
small-number statistics and are still at or beyond the cut- 
ting edge of current X-ray instruments. Different, com- 
plementary techniques must be employed to verify the 
phase membership of the O VI systems and confirm the 
presence of X-ray WHIM absorbers. 

The most convincing X-ray WHIM detections have 
relied on targeted searches of nearby sup ercluster fil- 
aments like the Sculp tor Wall detection (jBuote et al.l 
120091: IFang et all 12010ft . t he Pisces-Cetus Superc luster, 
and the Far Sculptor Wall (Zappaco sta et al.l l2010). The 
Sculptor Wall (at z = 0.030 ± 0.002 in the direction 
of the blazar target H 2356— 309) used both Chandra 
ACIS and XMM-Newton RGS to detect the O VII Ka 
absorption lin e at z = 0-032. The O VIII detection 
at z = 0.055 (IFang et al.ll2"00l towards PKS 2155-304 
is coincident with a group of narrow Lya lines at the 
redshift of a group of spiral galax ies (jShull et aJJ 11998b 
IShull. Timlinson. fc Girouxl [20031) . In many ways, the 
search for WHIM along the Mrk421 sight line (whether 
via X-ray or UV metal ions or BLAs) should also be con- 
sidered a targeted search, since this object lies just be- 
hind a section of the Great Wall, whether or not WHIM 
gas is ultimately confirmed. 

Although the current dataset appears to rule out the 
WHIM detections claimed by N05, due to the imprtance 
of these claims, we plan to re-observe Mrk421 with 
COS in HST Cycle 19. Deeper observations will ei- 
ther reveal a BLA feature below the current detection 
limit or place stronger constraints on WHIM gas in the 
sight line. Either way, the detection/non-detection of 
WHIM in the Great Wall will add an important data- 
point to the correlation of hot baryons and large-scale 
structure. The covering factor of WHIM gas around 
galaxy filamentary structure has significant impact on 
our understanding of the missing baryon problem at 
low-z. Correlatin g; the much larger catal ogs of low- 
z O VI detections (Danforth. fc Shulll 12008b iTripp et al.l 
20081 iThom fc Chenll2008D with galaxy redshift surveys 
may reveal additional filament- WHIM correspondence, 
and document the WHIM non-detections in similar LSS 
features. 
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